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ABSTRACT 

The  modeling  and  active  control  of  acoustic  pressure  in  a  2-D  cavity  with  a  flexible  bound¬ 
ary  (a  beam)  is  considered.  Control  is  implemented  in  the  model  via  piezoceramic  patches 
on  the  beam  which  are  excited  in  a  manner  so  as  to  produce  pure  bending  moments.  The 
incorporation  of  the  feedback  control  in  this  manner  leads  to  a  system  with  an  unbounded 
input  term.  Approximation  techniques  are  discussed  and  by  writing  the  resulting  system  as 
an  abstract  Cauchy  equation,  the  problem  of  reducing  interior  pressure  fluctuations  can  be 
posed  in  the  context  of  an  LQR  time  domain  state  space  formulation.  Examples  illustrating 
the  dynamic  behavior  of  the  coupled  system  as  well  as  demonstrating  the  viability  of  the 
control  method  on  a  variety  of  problems  with  periodic  forcing  functions  are  presented. 
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1  Introduction 


An  important  problem  in  the  study  of  structural  acoustics  involves  the  control  of  noise 
which  is  being  generated  by  the  vibration  of  an  elastic  structure.  An  example  of  this  is 
the  transmission  of  sound  through  an  airplane  fuselage  due  to  low  frequency  high  amplitude 
exterior  acoustic  fields.  The  need  for  lighter  materials  and  the  development  of  new  fuel 
efficient  turboprop  engines  has  reduced  the  transmission  loss  and  increased  the  acoustic  and 
vibration  levels  thus  causing  the  noise  control  problem  to  beconie  a  significant  impediment 
to  efficiency  In  recent  years. 

One  means  of  attempting  to  decrease  the  noise  levels  has  been  through  the  use  of  pas¬ 
sive  noise  control  techniques  which  n  ’  od  on  stiffened  structures  and  material  and  acoustic 
damping.  In  many  cases,  the  increa  d  weight  offset  the  advantages  obtained  through  the 
use  of  the  new  engines  and  lighter  materials. 

The  active  control  of  noise  in  this  setting  has  been  studied  both  in  a  frequency  domain 
setting  [11,  12,  13,  16,  17,  19]  and  from  an  infinite  dimensional  state  space  time  domain 
approach  (PDE  approach)  [1,  4,  5,  10].  Control  techniques  in  these  works  range  from  the  use 
of  acoustic  point  sources  which  generate  appropriate  secondary  pressure  waves  that  optimally 
interfere  with  the  offending  primary  pressure  wave  to  the  alteration  of  structural  dynamics 
through  applied  point  forces  and  bending  moments.  In  this  latter  case,  the  actuators  range 
from  electromagnetic  drivers  to  piezoceramic  patches. 

In  this  work,  we  continue  the  study  of  a  time  domain  state  space  formulation  in  which  the 
active  control  is  implemented  via  piezoceramic  patches  which  are  bonded  to  the  boundary 
of  the  acoustic  cavity  (a  more  theoretical  background  for  some  of  the  results  in  this  work 
can  be  found  in  [1]).  In  this  way,  we  can  take  advantage  of  the  natural  “feedback”  loop 
which  is  due  to  the  coupling  of  the  structural  vibrations  and  the  acoustic  fields.  Moreover, 
the  piezoceramic  patches  are  light,  can  be  bonded  directly  to  the  elastic  surface  and  are 
relatively  cheap  to  manufacture. 

The  example  we  consider  consists  of  an  exterior  noise  source  which  is  separated  from 
an  interior  chamber  by  an  elastic  plate.  This  plate  transmits  noise  or  vibrations  from  the 
exterior  field  to  the  interior  cavity  via  fluid/structure  interactions  thus  leading  to  the  for¬ 
mulation  of  a  system  of  partial  differentia!  equations  consisting  of  an  acoustic  wave  equation 
coupled  with  elasticity  equations  for  the  plate.  Control  is  implemented  in  the  example  via 
piezoceramic  patches  on  the  plate  which  are  excited  in  a  manner  so  as  to  produce  pure  bend¬ 
ing  moments  which  then  affect  the  bending  components  of  the  elasticity  equations.  Because 
the  incorporation  of  the  feedback  control  is  through  actuators  covering  only  sections  of  the 
boundary,  the  resulting  system  contains  an  unbounded  input  term.  Experiments  are  being 
designed  and  carried  out  at  ,N.\S,‘\  Langley  Research  Center  in  which  the  interior  cavity  is 
taken  to  be  cylindrical  with  a  circular  elastic  plate  a.iv!  sectorial  patches. 
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As  a  first  step  toward  developing  an  effective  linear  quadratic  regulator  (LQR)  state  space 
control  methodology  for  near  field  acoustic  problems  of  this  type,  it  is  useful  to  consider  a 
simplified  but  typical  model  consisting  of  a  2-D  interior  cavity  with  an  elastic  beam  at  one 
end  (see  Figure  1).  Here  T  represents  a  perturbing  force  on  the  beam  due  to  an  exterior 
noise  source.  This  in  turn  causes  fluctuations  in  the  interior  acoustic  pressure  field  and  hence 
unwanted  noise.  The  goal  in  the  control  problem  is  to  optimally  reduce  the  interior  pressure 
deviations  by  effecting  a  force  distribution  on  the  beam  that  decouples  the  cavity  acoustic 
response  from  the  beam  response  or  primary  excitation. 

In  Section  2,  a  model  set  of  differential  equations  for  the  problem  is  given  and  the  math¬ 
ematical  framework  needed  to  pose  the  control  system  in  variational  form  is  discussed.  The 
finite  dimensional  approximation  of  the  control  problem  is  presented  in  Section  3.  After 
discussing  suitable  choices  for  the  bases,  the  discrete  system  is  formulated  as  a  finite  dimen¬ 
sional  Cauchy  equation.  With  the  problem  in  this  format,  the  control  gains  can  be  obtained 
in  terms  of  the  solution  of  the  algebraic  Riccati  equation.  Section  4  contains  examples 
demonstrating  both  modeling  and  control  results.  The  first  example  illustrates  the  dynamic 
behavior  of  the  coupled  system  with  the  next  three  examples  demonstrating  the  viability  of 
the  control  method  on  a  variety  of  problems  with  periodic  forcing  functions. 


Figure  1.  Acoustic  chamber  with  piezoceramic  patches. 
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2  Mathematical  Model 


When  describing  the  interaction  between  the  acoustic  wave  motion  in  the  fluid  and  the 
vibration  of  the  structure,  both  the  fluid  velocity,  v,  and  the  acoustic  pressure,  p,  (the 
deviation  from  the  mean  pressure  at  equilibrium)  are  physically  significant.  In  order  to 
simplify  the  resulting  equations,  it  is  useful  to  introduce  a  velocity  potential  (p  which  is  a 
complex-valued  function  from  which  the  characteristic  values  v  and  p  can  be  calculated 
by  means  of  simple  differential  equations.  The  velocity  potential  is  determined  by  the 
relation  v{t,x,y)  =  ~'V(i){t,x,y)  and  the  pressure  is  related  to  this  velocity  potential  by 
p{t,x,y)  —  pf<f>t{t,x,y)  where  pj  is  the  equilibrium  density  of  the  fluid. 

For  acoustic  waves  with  small  amplitude,  both  the  potential  and  the  pressure  satisfy  the 
undamped  wave  equation  with  uniform  speed  of  sound  c  in  the  fluid  [14,  15];  hence 

=  c^A(j)  ix,y)  e  (l{t)  ,t  >  0  . 

The  boundaries  on  three  sides  of  the  variable  cavity  il[t)  are  taken  to  be  “hard”  walls  thus 
leading  to  the  zero  normal  velocity  boundary  conditions 

V^-n  =  0  (x,j/)€r,t>0 

where  n  is  the  outer  normal.  It  is  assumed  that  the  perturbable  boundary  consists  of  an 
impenetrable  fixed-end  Euler-Bernoulli  beam  with  Kelvin- Voigt  damping.  If  w{t,  x)  is  used 
to  denote  the  transverse  displacement  of  the  beam  with  linear  mass  density  pt,  the  equations 
of  motion  are 


PbWtt  +  =  -pf<l)t{t,x,u'{t,x))  -I-  f{t,x) 

^XjU  div 

wit,0)  = =  w{t,a)  = =  0  i>0  , 


0  <  X  <  a  , 
f  >  0  , 


where  M{t,  x)  is  the  internal  moment  and  /  is  the  external  applied  force  due  to  pressure  from 
an  exterior  noise  field.  For  an  uncontrolled  beam  with  Kelvin- Voigt  damping,  the  moment 
contains  both  strain  and  strain  rate  components  and  is  given  by 

«  f/i  \  f-. 

Mit,x)-  EI~  +  cdI^^^  . 

The  final  coupling  equation  is  the  continuity  of  velocity  condition 


Wt(t,x)  =  V4>(t,x,w(t,x)}  •  n  ,  0  <  X  <  a  >  0 


(2.2) 


which  results  from  the  assumption  that  the  beam  is  impenetrable  to  fluid.  Under  an  as¬ 
sumption  v/f  small  displacements  (w(t,x)  =  w(t,x)  -|-  S  where  =  0)  which  is  inherent  in 
the  Euler-Bernoulli  formulation,  the  beam  equation  in  (2.1)  can  be  approxim.*tcd  by 


Pb^tt  "f 


g2 

—  M(t,x)  =  -Pf[<f>t(t,x,0)  +  -h  /(f,x) 
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while  (2.2)  can  be  approximated  by 


Wt{t,  x)  —  a;,0)  •  h  +  (V^y(<,  x,  0)re)  •  h  . 


To  first  order,  these  last  two  equations  can  be  approximated  by  dropping  the  higher  order 
terms  —pj<f>ty{t,x,0)w  and  {V(l)y{t,x,0)w)  ■  h.  Then  upon  approximating  the  domain  n(t) 
by  the  fixed  domain  Q  =  [0,a]  x  [0,^],  we  obtain  the  approximate  uncontrolled  model 


A(?  (r,  y)  e  fl  ,  f  >  0  , 

^(p-n  —  O  {x,  y)  G  r  ,  t  >  0  , 

=  — u’((/,x)  0  <  a:  <  «  >  0  , 

dy 

(  d^w  d^w  N 

+  fit, X) 

w{t,^)  =  -  w{t^a)  -  —{t,a)  -  Q  t>0  , 


(2.3) 


(?(0,x,y)  =  (;5o(x,y)  ,  tu(0,x)  =  u;o(x) 

(?l)((0,x,y)  =  <pi(x,y)  ,  u;((0,x)  =  u;i(x)  . 


For  control  of  structural  vibrations  and  the  acoustic  pressure  field  in  this  model,  s  piezo¬ 
ceramic  patches  are  attached  to  the  beam  as  shown  in  Figure  1.  These  patches  are  excited 
in  a  manner  so  as  to  produce  pure  bending  moments  ([6,  7,  9])  (see  Figure  2).  If  H  is  used 
to  denote  the  Heaviside  function,  the  model  for  the  controlled  beam  can  be  written  as 


^ 


+  CpI 


d^w  '\ 
dx'^dt ) 


T  P j 1  ^ 


d^ 

dx^ 


T 


^u,(f)  [H{x  -  a,i)  -  H{x  -  0,2)] 


+  f{t,x)  . 


(2.4) 


Here  u,{t)  is  the  voltage  applied  to  the  patch,  is  a  parameter  which  depends  on 
the  geometry  and  piezoceramic  material  properties,  T  is  the  patch  thickness  and  d^i  is  the 
piezoelectric  strain  constant  (sec  [6,  7]).  It  should  be  noted  that  the  incorporation  of  (2.4) 
into  (2.3)  leads  to  a  system  with  an  unbounded  input  term  since  it  involves  the  second 
derivative  of  the  Heaviside  function. 
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Figure  2.  Piezoceramic  patch  excitation. 


In  order  to  pose  this  problem  in  a  manner  which  is  conducive  to  approximation  and 
control,  the  state  is  taken  to  be  =  {(/>■,  w)  in  the  state  space  IJ  =  X  Here 

is  the  quotient  space  of  over  the  constant  functions.  The  use  of  the  quotient  space 
results  from  the  fact  that  the  potentials  are  determined  only  up  to  a  constant. 

To  provide  a  class  of  functions  which  are  considered  when  defining  a  variational  form 
of  the  problem,  we  also  define  the  Hilbert  space  V  =  x  Hq{To)  where  is 

the  quotient  space  of  IP  over  the  constant  functions  and  //^(Tq)  is  given  by  Hq{To)  = 
{0  €  lP{To)  :  0(a')  =  'p{x)  =  0  at  x  =  0,a}. 

A  complete  discussion  concerning  the  first  and  second  order  weak  or  variational  forms  of 
the  control  system  is  given  in  [1].  For  our  purposes  here,  it  suffices  to  note  that  from  (2.3) 
and  (2.4),  it  follows  that  the  second-order  system  in  variational  form  is  given  by 


+  /  PbWttrid'y 

JVo 

+  /  -f  f  El D^wD^nd'y 

-f  f  {cdI D'^WtD'^r]  +  pj{(i>tr]  -  Wt^)}d-y 

'•'Vq 


/K  ^  (I  ^ 

0  ■*1=1 


(2.5) 


+ 


for  all  (^,f/)  in  V.  In  order  to  simplify  the  above  expression,  we  have  adapted  the  notation 
//,j(x)  =  H{x  —  Utj) ,  i  =  1,  ■  ■  ■  ,s,  y  =  1,2.  The  system  can  be  written  in  first-order  form 
by  defining  the  product  spaces  V  =  V  x  V  and  Ti,  =  V  x  H  and  taking  the  state  to  be 
Z  —  (0,  uj,  0,  tk] .  Note  that  the  state  now  contains  a  multiple  of  the  pressure  since  p  —  pj^. 
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The  first-order  variational  form  is 


/  +  /  Pb{w)tTjd'y 

C  JVq 


+ 

+ 


/  pjV<f>  ■  Vidu;  +  /  EID^wD^rjd-f 
Jq  Jvq 

f  {cdID'^wD'^t]  +  pf{^r]  -  w^)}dj 

^To 

1 


+ 


(2.6) 


for  all  [iiTi)  in  V.  Again,  a  more  complete  discussion  concerning  the  formulation  of  the 
first-order  system  in  weak  form  is  given  in  [1]. 


3  Finite  Dimensional  Approximation  and  Control 

The  problem  (2.6)  is  now  in  a  form  which  is  amenable  both  to  approximation  and  control.  As 
shown  in  [1],  the  infinite  dimensional  system  can  be  formally  written  as  an  abstract  Cauchy 
equation,  and  periodic  infinite  dimensional  control  results  similar  to  those  found  in  [8]  can 
be  applied  (see  also  [2,  3]).  Approximate  feedback  gains  (and  hence  controlling  voltages)  can 
be  obtained  by  using  a  standard  Galerkin  approach  in  which  one  chooses  a  sequence  of  finite 
dimensional  subspaces  C  H  with  projections  :  Ti  — >  .  In  order  to  guarantee 
the  convergence  of  the  approximate  gains  to  those  of  the  infinite  dimensional  system,  it  is 
sufficient  to  impose  various  conditions  on  the  original  and  approximation  systems.  These 
hypotheses  include  convergence  requirements  for  the  uncontrolled  problem  as  well  as  the  re¬ 
quirement  that  the  approximation  systems  preserve  stabilizability  and  detectability  margins 
uniformly. 

Cubic  splines  are  used  as  a  basis  for  the  beam  since  they  satisfy  the  smoothness  re¬ 
quirements  as  well  as  being  easily  implemented  when  adapting  to  the  fixed-end  boundary 
conditions  and  patch  discretizations.  Letting  denote  the  cubic  splines  which  have 

been  modified  to  satisfy  the  boundary  conditions  (see  [1]  for  details),  the  corresponding  n  —  1 
dimensional  approximating  subspace  is  given  by  //"  =  span{5"}"rj^  and  the  approximate 
beam  solution  is  taken  to  be 


»=i 

The  2-D  cavity  basis  is  taken  to  be  ,  where  B"‘(x,j/)  is  the  tensor  product  of  1-D 

Legendre  polynomials  which  have  been  scaled  by  transformation  to  the  x  and  y  intervals  [0,  a] 
and  [0,  £],  respectively.  The  total  number  of  cavity  basis  functions,  m  =  (mx  +  l)-(my-|-l)  — 1, 
reflects  the  deletion  of  the  constant  function  so  as  to  guarantee  that  the  set  of  functions  is 
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suitable  as  a  basis  for  the  quotient  space.  No  boundary  modifications  are  needed  for  the 
Legendre  basis  since  natural  boundary  conditions  occur  on  all  sides  of  the  cavity.  The 
m  dimensional  approximating  subspace  is  then  taken  to  be  //"*  =  span{ii™}^j  and  the 
approximate  cavity  solution  is  given  by 

m 

i=l 

The  approximating  state  space  is  then  taken  to  be  x  where  N  =  m  +  n  —  1, 

and  the  product  space  for  the  first  order  system  is  x  .  As  shown  in  [1],  the 

restriction  of  the  infinite  dimensional  system  (2.5)  to  7i^  x  yields 

=  A^y^{t)  +  B^u{t)  +  F^{t) 


or  equivalently  the  finite  dimensional  Cauchy  equation 


y‘'^it)  =  A^y^it)  +  +  F‘^{t) 

j/''^(0)  =  y^  . 


Explicit  descriptions  of  the  mass  and  stiffness  matrices  AF'^  and  as  well  as  detailed  defi¬ 
nitions  of  the  control  matrix  B"  cind  the  force  vector  F^{t)  can  be  found  in  [1].  The  2N  x  1 

vector  y^{t)  =  (<(0, '  •  • ,  •  •  • ,  •  •  • ,  {t) ,  •  •  • , 

contains  the  approximate  state  coefficients  while  u(t)  =  (ui(t),  •  •  •  contains  the  s  con¬ 

trol  variables.  The  periodic  finite  dimensional  control  problem  is  then  to  find  u  6  L^(0,t) 
which  minimizes 


=  \l^  {(gV(0,2/"'(0)F-  +  {Ru{t),u{t))f,s}dt  ,  iV  =  m  -f  n  -  1 

where  y^  solves  (3.1),  r  is  the  period,  /?  is  an  s  x  s  diagonal  matrix  and  r„  >  0,  f  =  1,  •  •  • ,  s 
is  the  weight  on  the  controlling  voltage  into  the  patch. 

The  nonnegative  definite  matrix  can  be  chosen  so  as  to  emphasize  the  minimization 
of  particular  state  variables  as  well  as  to  create  windows  that  can  be  used  to  decrease 
state  variations  of  certain  frequencies.  From  energy  considerations  as  discussed  in  [1],  an 
appropriate  choice  for  in  this  case  is 


=  M^V 


where  again  is  the  mass  matrix,  and  the  diagonal  matrix  T)  is  given  by 

T>  =  diag  \d,r\d2r-\d^r,d^r-^^ 


Here  F  ,  k  =  m,n  —  \  ,  denotes  a.  k  x  k  identity  and  the  parameters  d,  are  chosen  to  enhance 
stability  and  performance  of  the  feedback. 

The  optimal  control  is  then  given  by 


( 


u^{t)  =  R-\B^f 


where  FI'''  is  the  solution  to  the  algebraic  Riccati  equation 


)rii.v  +  Q‘'''  =  0  . 


(3.2) 


F\ir  the  regulator  problem  with  periodic  forcing  function  r^{t)  solves  the  linear  dif¬ 

ferential  equation 


r^'(O)  =  r^\r) 


(3.3) 


while  the  optimal  trajectory  is  the  solution  to  the  linear  differential  equation 

y'"’(0)  =  -/'^'(r)  • 


4  Numerical  Results 


The  general  problem  under  consideration  in  the  following  examples  is 

p((  =  c^Ap  (.i'lf/)  €  n  ,f  >  0  , 

vp •  n  =  u  (j-, (/)  e  r  ,<  >  0 , 

(Jo 

(Ti',0)  =  —LCt{t.x)  0  <  X  <  .6  ,  <  >  0  , 

oy 


Pb^'tt  + 


0^ 


dx^  dxm 


>)2  /  /v^(7  \ 

=  ^  (^^'^~7-'~^(0[^^(-^  -  »ii)  -  -  Qi'2)]j 


-  pf(f)t{t,x,Q)  +  f{t,x) 


Lv{t,0)  =  ^(<,0)  =  w{t,.6)  -  ^(f,.6)  =  0 
ax  ax 


0  <  X  <  .6  ,  f  >  0  , 
f  >  0  , 


(4.1) 


p(0,x,y)  =  p,(0,x  y)  =  u;(0,x)  =  n;t(0,x)  =  0  . 

'Fhe  parameter  choices  a  —  .()  ui,  f  =  1  m,  pj  =  1.21  kg/nF,  cJ  =  117649  m'^(sec^, 
Pb  =  \.'J5k(j/rn,  FI  =  73.96  XnF,  cpI  =  .001  kgriF/sec,  A'®  =  82.9629,  T  =  .0005  m, 
and  dsi  =  1.9  x  10“‘°  ai/V ,  are  physically  reasonable  for  a  .6  rn  by  1  rn  cavity  in  which  the 
bounding  end  beam  has  a  centered  piezoceramic  patch.  The  beam  is  assumed  to  have  width 
and  thickness  .1  rn  and  .005  ni,  respectively.  For  an  undamped  fixed-fixed  beam  with  the 
above  dimensions  and  density,  the  natural  frequencies  of  the  first  three  modes  are 

/i  =  73.2  hertz  , 

Ji  —  201.8  hertz  , 

/a  =  395.6  hertz  , 
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respectively.  Likewise,  for  pressure  oscillations  in  a  cavity  of  this  size,  the  natural  frequencies 
of  the  modes  are  given  by  ^  where  L  =  C  =  I  in  the  y-direction  and 

Z-  =  a  =  .6  in  the  x-direction.  .Note  that  Loth  sets  of  frequencies  are  for  the  uncoupled 
systems. 


4.1  System  Dynamics 

In  this  example,  we  deicrmine  the  behavior  of  the  uncontrolled  coupled  system  when  the 
forcing  function  /  is  taken  to  be  an  impulse  at  the  center  of  the  beam;  that  is, 

/(t,x)  =  6{x  -  .3]S{t)  . 

F’roin  the  response,  the  natural  freciuencies  and  dynamic  behavior  of  the  system  can  be 
determined  and  compared  to  those  of  the  uncouuled  beam  and  cavity.  .Note  that  from  the 
definitions  of  the  component  matrices  and  vectors,  it  follows  immediately  that  the  right  hand 
side  vector  has  elements  given  by 

[efiol  =  / 

J\  0 

where  denotes  a  modified  cubic  spline  basis  function.  Thus  F2  {t)  is  simply  given  by 
^'2  i^)  —  (0)  ■  ■  •  )b,  1,-1,  1,0,  •  •  •  ,0]^  S{t).  .Note  that  this  contribution  is  included  only  at  lime 
T  =  0  after  which  the  system  is  allowed  to  run  unforced  through  time  T  =  20115.  This 
temporal  interval  is  sufiicienl  for  demonstrating  the  dynamic  behavior  of  the  system. 

The  beam  and  pressure  responses  obtained  with  =  lUy  —  12  and  n  =  IG  basis  func¬ 
tions  at  the  points  .V  =  .0  and  (-V,  y)  =  (.3,  .1)  are  plotted  in  Figure  3  with  corresponding 
frequency  plots  in  Figure  -1.  From  Figure  3  it  can  be  seen  that  due  to  the  damping  in  the 
beam,  both  the  beam  and  pressure  oscillations  are  gradually  decreasing  in  magnitude  as  the 
system  is  allowed  to  run  unforced.  I'he  frequency  plots  in  Figure  4  indicate  the  presence 
of  several  natural  frequencies  for  the  coupled  system.  Although  one  would  expect  these 
frequencies  to  correspond  in  some  manner  to  the  individual  natural  frequencies  of  the  beam 
and  cavity,  one  must  also  note  that  the  system  (4.1)  involves  not  only  coupling  between 
the  cavity  and  beam  equations  but  also  includes  damping  in  the  beam.  Indeed,  the  system 
responses  at  65.9  and  387.8  hertz  are  slightly  lower  in  frequency  than  the  natural  frequencies 
of  the  first  two  symmetric  modes  of  the  undamped  beam  which  have  values  of  73.2  and  395.6 
hertz.  The  remaining  resonant  peaks  are  due  to  pressure  oscillations  in  the  cavity.  The 
responses  at  181.3  and  343.9  hertz  arc  slightly  higher  than  the  values  of  171.5  and  343  :  rtz 
which  are  the  first  two  natural  frequencies  in  the  y-direction  for  pressure  oscillations  in  a 
cavity  of  this  size.  .A  response  at  519  hertz  corresponding  to  the  third  harmonic  (514  hertz) 
in  the  y-direction  ca.i  be  seen  if  one  considers  data  taken  over  the  much  shorter  time  inter¬ 
val,  [0,3/75].  This  response  has  very  low  energy  however,  and  has  died  away  long  before 
T  =  20/75.  As  seen  in  Figures  11  and  18  of  Sections  4.3  and  4.4,  system  oscillations  at  this 
frequency  are  more  pronounced  when  produced  by  a  forcing  function  w'hich  is  acting  over  a 
longer  temporal  interval. 
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The  system  response  at  572.4  hertz  corresponds  to  the  uncoupled  cavity  frequency  of 
571.7  hertz  for  the  second  x  mode,  f^ecause  the  excitation  is  symmetric  in  x,  one  does  not 
see  a  response  corresponding  to  the  nonsyrnmetric  first  x  mode  having  the  frequency  285.8 
hertz.  It  is  more  difficult  to  determine  correspondences  between  higher  frequency  system 
modes  (600  hertz  and  abo\’e)  and  those  of  the  uncoupled  beam  and  cavity.  This  may  in 
part  be  due  to  higher  frequency  effects  of  coupling  and  damping.  Finally,  it  should  be  noted 
that  the  mode  shapes  for  the  coupled  system  as  well  as  the  natural  frequencies  are  slightly 
different  than  those  of  the  component  beam  and  cavity;  hence  care  must  be  taken  when 
describing  the  dynamics  of  the  coupled  system  in  terms  of  the  properties  of  the  undamped 
and  uncoupled  beam  and  cavity. 


X 10'^^  L’ncon(roUe<l  Heani  Dispi«cemen(  at  X  =  .3 


Figure  3.  The  beam  and  pressure  responses  to  a  centered  impact. 
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Figure  4.  I'he  beam  and  pressure  frequency  responses  to  a  centered  impact. 
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4.2  Resonant  Excitation 


In  the  following  examples,  tne  effects  of  feedback  control  on  the  system  (4.1)  for  various 
forcing  functions  are  described.  The  forcing  function  in  this  example  is 

f{t,x)  —  2.04sin(1507r<) 

which  corresponds  to  a  frequency  of  75  hertz  and  hence  is  close  to  the  natural  frequency  of 
the  first  system  mode  which  is  65.9  hertz  (see  Section  4.1).  This  models  a  periodic  exterior 
plane  wave  with  a  root  mean  square  (rms)  sound  pressure  level  of  117  dB. 

Control  was  implemented  via  a  single  centered  piezoceramic  patch  covering  1/6  of  the 
beam  length  (see  Figure  5).  The  quadratic  cost  functional  parameters  were  taken  to  be 
di  =  ^2  =  ^4  =  1,  <^3  =  10“*  and  R  =  10“®  with  of  much  larger  magnitude  than  di,d2 
or  ^4  to  emphasize  the  penalization  of  large  pressure  variations.  When  running  the  uncon¬ 
trolled  problem,  the  beam  and  pressure  oscillations  could  be  fully  resolved  with  the  choices 
rrij.  =  rriy  =  4  and  n  —  8.  In  the  control  problem  however,  a  larger  number  of  basis  functions 
was  needed  to  fully  resolve  the  early  transient  behavior  and  the  results  which  follow  were 
obtained  with  the  choices  nij.  =  niy  =  8  and  n  —  12.  Note  that  this  results  in  a  total  of  80 
cavity  and  11  beam  basis  functions. 


0  .6 


Figure  5.  Acoustic  chamber  with  one  centered  1/6  length  piezoceramic  patch. 

In  order  to  solve  for  the  optimal  control  and  trajectory,  it  is  necessary  to  solve  not  only 
the  algebraic  Riccati  equation  (3.2),  but  also  the  trajectory  equation  (3.4)  and  the  tracking 
equation  (3.3).  The  gains  from  the  Riccati  equation  were  calculated  via  Potter’s  method 
(see  [18]).  Because  numerical  evidence  indicated  that  the  unconstrained  solutions  to  the 
trajectory  and  tracking  equations  were  roughly  periodic  with  period  r  =  1/75,  the  problems 
were  solved  as  initial  value  problems  with  starting  values  y{0)  =  0  and  r(  10/75)  =  0  rather 
than  as  free  boundary  value  problems.  The  choice  for  initial  state  is  physically  reasonable 


12 


while  the  choice  to  integrate  backwards  in  time  in  (3.3)  is  made  to  reduce  numerical  instability 
when  solving  the  ODE  system  for  r^{t). 

The  uncontrolled  and  controlled  approximate  beam  displacements  and  acoustic  pressures 
{p^  =  pj<p^)  at  the  points  A”  =  .3  and  {X,Y)  =  (.3,  .1)  are  plotted  in  Figure  6  for  the 
time  interval  [0,10/75].  The  uncontrolled  solutions  exhibit  a  beat  phenomenon  which  re¬ 
sults  from  the  fact  that  the  frequency  of  the  forcing  function  is  slightly  greater  than  the 
natural  frequency  of  the  first  system  mode.  From  this  data,  it  can  be  determined  that  the 
uncontrolled  rms  pressure  at  the  point  {X,Y)  =  (.3,  .1)  is  pma  =  1.0957  N/m"^.  This  yields 
an  uncontrolled  interior  sound  pressure  level  of  94.6  dB.  After  a  transient  interval,  the  con¬ 
trolled  solutions  are  periodic  and  are  maintained  at  a  level  which  is  approximately  10%  of 
that  found  in  the  uncontrolled  case.  The  controlled  rms  pressure  is  prms  =  -125  N/m'^  and 
the  controlled  interior  sound  pressure  level  is  maintained  at  75.8  dB  which  is  a  18.8  dB 
reduction. 

The  frequencies  of  the  uncontrolled  and  controlled  beam  displacements  and  cavity  pres¬ 
sures  are  plotted  in  Figure  7  (note  the  differing  scales  in  the  plots).  The  uncontrolled  beam 
plot  shows  the  presence  of  the  forcing  frequency  as  well  as  the  slightly  lower  natural  system 
frequency  while  the  uncontrolled  cavity  plot  also  shows  a  slight  system  response  at  180  hertz. 
The  controlled  response  in  both  cases  shows  only  the  presence  of  the  forcing  frequency. 

To  further  illustrate  the  state  reduction  with  feedback  control,  the  uncontrolled  and 
controlled  acoustic  pressures  at  the  times  T  =  1/75,  2/75,  6/75  and  10/75  are  plotted  in 
Figure  9.  These  are  shown  as  surface  plots  illustrating  the  pressure  amplitude  distribution 
throughout  the  acoustic  cavity.  These  results  are  representative  of  those  found  throughout 
the  time  interval  [0, 10/75]  and  in  conjunction  with  Figure  6,  demonstrate  that  the  pres¬ 
sure  and  beam  displacement  are  uniformly  reduced  and  maintained  at  a  very  low  level  of 
magnitude  in  spite  of  the  periodic  forcing  function. 

The  controlling  voltage  u(t)  is  plotted  in  Figure  8.  As  expected,  it  is  periodic  with 
period  1/75.  It  should  be  noted  that  the  magnitude  of  u{t)  remains  less  than  60F  which  is 
a  physically  reasonable  voltage  to  put  into  the  piezoceramic  patches. 


13 


(NAn^2)  DispUcement  (m) 


0  0.02  0.04  0.06  0.08  0.1  0.12  0  0.02  0.04  0.06  0.08  0.1  0.12 

Time  (sec)  Time  (sec) 


Time  (sec)  Time  (sec) 


Figure  6.  Uncontrolled  and  controlled  beam  displacements  and  pressures  at  the  points 
X  =  .3  and  (X,  Y)  =  (.3,  .1)  throughout  the  time  interval  [0, 10/75]. 
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CoQirolkd  PreMuie 


4.3  Multiple  Resonant  Excitation 

The  forcing  function  for  this  example  is 

f{t,  x)  —  2.04  [sin(1507rt)  +  sin(7907rt)] 

and  hence  excites  the  first  and  fourth  system  modes  which  have  natural  frequencies  of  65.9 
and  387.8  hertz,  respectively  (see  Section  4.1).  This  models  a  periodic  plane  wave  with  a 
sound  pressure  level  of  120  dB. 

Control  was  again  implemented  via  a  single  centered  piezoceramic  patch  covering  1/6 
of  the  beam  length  as  shown  in  Figure  5.  The  results  reported  below  were  obtained  with 
the  choices  di  =  d2  =  =  I,  =  10^  and  R  =  10“*^  for  the  quadratic  cost  functional 

parameters.  Several  other  combinations  were  tried  but  these  reported  values  provided  a  good 
balance  between  the  amount  of  pressure  reduction  obtained  and  the  amount  of  controlling 
voltage  applied  to  the  piezoceramic  patches. 

From  numerical  tests,  it  was  found  that  the  uncontrolled  and  controlled  beam  and  pres¬ 
sure  oscillations  could  be  resolved  with  the  choices  rUx  =  niy  =  10  and  n  =  16  and  the 
following  results  were  obtained  with  these  values.  The  need  for  a  larger  number  of  basis 
functions  than  were  used  in  Section  4.2  reflects  the  presence  of  higher  frequency  oscillations 
which  resulted  from  the  higher  frequency  component  of  the  forcing  function. 

The  time  interval  of  interest  in  this  example  was  taken  to  be  [0, 9/75]  and  the  uncontrolled 
beam  displacements  and  acoustic  pressures  at  the  points  X  =  .3  and  {X^Y)  =  (.3,.l)  are 
plotted  in  Figures  10  and  14.  Control  was  implemented  at  time  T  =  0  to  obtain  the 
controlled  solution  plots  in  Figure  10  while  the  corresponding  plots  in  Figure  14  show  the 
results  obtained  when  the  system  was  allowed  to  run  uncontrolled  until  T  =  3/75  =  .04 
seconds  at  which  point  the  controlling  voltage  was  applied.  The  frequency  plots  in  Figure  11 
correspond  to  the  results  in  Figure  10.  The  primary  frequencies  in  both  the  uncontrolled 
and  controlled  problems  are  the  driving  frequencies  of  75  and  395  hertz.  In  the  uncontrolled 
case,  there  is  also  a  low  energy  system  response  at  519.5  hertz  (recall  that  the  third  harmonic 
in  the  j/-direction  for  the  uncoupled  cavity  is  at  514  hertz).  By  comparing  the  scales  in  the 
plots  of  Figure  11,  one  can  see  that  the  energy  in  the  controlled  frequency  responses  is 
approximately  one  tenth  of  that  found  in  the  uncontrolled  responses. 

From  Figure  10,  it  is  seen  that  when  control  is  started  at  T  =  0,  the  solutions  are 
periodic  and  are  maintained  at  a  level  which  is  about  12%  of  that  found  in  the  uncontrolled 
case.  By  calculating  the  root  mean  square  pressures,  it  was  determined  that  at  the  point 
(X,  F)  =  (.3,  .1),  the  uncontrolled  sound  pressure  level  is  97.75  dB  whereas  the  controlled 
sound  pressure  in  this  case  is  maintained  at  a  level  of  78.93  dB.  The  two  dimensional  surface 
plots  in  Figure  13  show  spatial  slices  of  the  uncontrolled  and  controlled  pressures  at  the 
times  T  =  1/75,3/75,5/74  and  9/75  for  the  case  when  control  was  started  at  T  =  0.  These 
plots  in  conjunction  with  those  in  Figure  10  show  that  even  in  the  case  of  strong  multiple 
frequency  input,  the  pressure  and  beam  displacement  are  uniformly  reduced  and  maintained 
at  a  very  low  level  throughout  the  time  period.  They  also  demonstrate  that  the  small  near 
field  acoustic  response  which  remains  near  the  beam  quickly  decays  as  one  moves  further 
into  the  cavity. 

The  controlling  voltage  u{t)  obtained  with  control  starting  at  T  =  0  is  plotted  in 
Figure  12.  As  in  Section  4.2,  it  is  periodic  although  in  this  case,  it  also  reflects  the  high 
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frequency  contributions  of  the  forcing  function.  The  magnitude  remains  below  70K  which 
again,  is  a  physically  reasonable  voltage  to  apply  into  the  piezoceramic  patch.  We  remark 
that  this  voltage  can  be  decreased  without  any  loss  of  control  by  using  longer  centered 
patches. 

As  seen  in  Figure  14,  the  implementation  of  control  at  T  =  3/75  =  .04  seconds  implies 
that  the  control  scheme  is  starting  with  an  initial  beam  displacement  of  —2.4  x  10“®  m  at 
A"  =  .3  and  an  initial  pressure  of  1.8  N/rn^  at  (A,  Y)  =  (.3,  .1).  By  T  =  .06  however,  both 
the  beam  displacement  and  pressure  have  been  reduced  to  the  levels  obtained  when  the 
controlling  voltage  was  applied  rd  T  =  0  (see  Figure  10),  and  the  solutions  are  maintained  at 
this  level  throughout  the  rest  of  the  time  interval.  The  controlling  voltage  u{t)  for  this  case 
is  plotted  in  Figure  15.  As  expected,  it  is  quite  large  when  the  control  is  initially  applied  but 
quickly  decreases  and  is  periodic  with  magnitude  less  than  701^  throughout  the  second  half 
of  the  time  interval.  The  magnitude  of  the  controlling  voltage  can  again  be  reduced  by  using 
longer  centered  patches.  These  results  demonstrate  that  this  methodology  is  useful  not  only 
for  maintaining  a  low  level  of  pressure  oscillations  throughout  the  time  interval,  but  also  for 
reducing  pressure  oscillations  which  have  been  allowed  to  build  up  in  the  cavity. 


Figure  10.  Uncontrolled  and  controlled  beam  displacements  and  pressures  at  the  points 
A  =  .3  and  (A,  Y)  =  (.3,  .1).  Control  is  started  at  T  =  0. 
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Figure  11.  Uncontrolled  and  controlled  beam  and  cavity  frequencies  when  control  is  started  at  T  =  0. 


The  Optimal  Control 


Figure  12.  The  optimal  control  u{t)  when  control  is  started  at  T  =  0. 
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Figure  15.  The  optimal  control  u(t)  when  control  is  started  at  T  =  3/75  =  .04  seconds. 


4.4  OfF-Resonant  Excitation 

In  this  example,  the  forcing  function  was  taken  to  be 

f{t,x)  =  2.04  sin(4707rt) 

with  a  frequency  of  235  liertz  which  is  half  way  between  the  natural  frequencies  of  the  first 
and  fourth  system  modes.  As  noted  in  Section  4.2,  this  models  a  periodic  plane  wave  with 
an  rms  sound  pressure  level  of  117  dB. 

Control  in  this  example  was  implemented  via  a  single  centered  piezoceramic  patch  cov¬ 
ering  1/2  of  the  beam  length  as  shown  in  Figure  16.  As  in  the  previous  examples,  sev¬ 
eral  combinations  of  the  quadratic  cost  functional  parameters  were  tried  with  the  choices 
di  =  ^2  =  d4  =  1,  da  =  lO'*  and  R  =  10“®  providing  a  good  balance  between  the  amount 
of  controlling  voltage  applied  to  the  piezoceramic  patches  and  the  amount  of  reduction  in 
pressure  fluctuations  which  could  be  obtained.  Multiple  patches  of  shorter  length  were  also 
tried  but  it  was  found  that  for  the  symmetric  forcing  function  being  applied,  the  single  patch 
provided  the  best  control. 


Figure  16.  Acoustic  chamber  with  one  centered  1/2  length  piezoceramic  patch. 

From  numerical  tests,  it  was  determined  that  the  uncontrolled  solutions  could  be  resolved 
with  rux  =  rriy  =  8,  n  =  12  basis  functions  whereas  rrix  =  =  12,  n  =  16  were  needed  to 

resolve  the  controlled  beam  and  pressure  oscillations  due  to  the  transient  high  frequency  re¬ 
sponses  to  the  controlling  voltage.  The  results  below  were  obtained  with  these  discretization 
sizes. 

From  Figure  18,  it  can  be  seen  that  the  uncontrolled  responses  exhibit  not  only  the 
driving  frequency  but  also  transient  excitations  at  65.9,  181.6,  345.2,  387.7  and  519.5  hertz 
which  are  due  to  the  natural  frequencies  of  the  coupled  system.  The  presence  of  the  multiple 
frequencies  can  also  be  seen  in  the  plots  of  the  uncontrolled  beam  displacement  and  pressure 
which  are  given  in  Figure  17.  The  time  interval  in  both  cases  was  taken  to  be  [0, 16/235]. 
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The  control  plots  in  Figure  17  demonstrate  that  after  a  transient  phose  of  about  three 
periods,  both  the  beam  displacement  and  pressure  are  significantly  reduced  and  maintained 
at  a  low  magnitude  throughout  the  rest  of  the  time  interval.  At  the  point  =  (-3, -1), 

the  uncontrolled  sound  pressure  level  is  82.8  dB  whereas  the  controlled  sound  pres-ure  level  is 
reduced  15.7  dB  to  67.1  dB.  The  frequency  plots  of  the  controlled  solutions  (see  Figure  18) 
show  that  the  dominant  response  is  the  driving  frequency  of  235  hertz.  These  plots  also 
indicate  the  presence  of  a  higher  frequency  response  to  the  control  at  519.5  hertz  which  is 
much  more  significant  than  in  the  uncontrolled  case.  This  response  is  transient  however,  and 
results  taken  further  out  in  time  indicate  that  by  approximately  48  periods,  ii  is  gone.  We 
have  also  found  that  when  longer  patches  are  used  (for  example,  2/3  length  patches),  the 
transient  high  frequency  response  is  much  smaller  than  that  observed  here. 

To  further  illustrate  the  state  reduction  with  feedback  control  via  the  piezocer''mic 
patches,  the  uncontrolled  and  controlled  acoustic  pressur<’s  at  the  times  T  =  1/235, 
T  =  3/235,  T  =  10/235  and  T  =  16/235  are  plotted  in  Figure  20.  These  results  reinforce 
those  in  shown  in  Figure  17  which  show  that  after  a  transient  phase  in  which  the  controlled 
pressure  undergoes  high  frequency,  low  magnitude  oscillations,  the  pressure  field  settles  into 
the  frequency  of  the  driving  force  and  is  uniformly  maintained  at  a  low  level  throughout  the 
time  interval. 

Figure  19  contains  a  plot  of  the  controlling  voltage,  u{t).  It  too  is  periodic  and  remains 
at  a  magnitude  of  less  than  25  volts. 


Figure  17.  Uncontrolled  and  controlled  beam  displacements  and  pressures  at  the  points 
A  =  .3  and  (X,  Y)  =  (.3,  .1)  throughout  the  time  interval  [0, 16/235]. 


23 


*10*^  Uncontrolled  Beam  Frequency  ai  X  -  .3 
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Figure  18.  Uricoalrollecl  and  controlled  beam  and  cavity  frequencies. 


The  Opiimal  Control 


Figure  19.  The  optimal  control  ti(0- 


Figure  20.  Uncontrolled  and  controlled  pressures. 


5  Conclusion 


In  this  work  we  continued  the  study  of  active  noise  control  techniques  for  2-D  structural 
acoustics  problems.  The  example  which  we  considered  consists  of  a  2-D  interior  cavity  with 
an  active  beam  at  one  end.  A  perturbing  force  due  to  an  exterior  noise  source  is  placed  on 
the  beam  thus  causing  fluctuations  in  the  interior  acoustic  pressure  field  and  hence  unwanted 
noise.  Control  is  implemented  in  the  model  via  piezoceramic  patches  on  the  beam  which  are 
excited  in  a  manner  so  as  to  produce  pure  bending  moments.  Approximation  techniques  are 
discussed  and  by  writing  the  resulting  system  as  an  abstract  Cauchy  equation,  the  problem 
of  reducing  interior  pressure  fluctuations  can  be  posed  in  the  context  of  an  LQR  time  domain 
state  space  formulation. 

Several  examples  are  considered  which  demonstrate  both  modeling  and  control  results. 
The  first  example  illustrates  the  dynamic  behavior  of  the  coupled  system.  By  looking  at 
the  natural  frequencies,  it  can  be  seen  that  although  many  of  the  features  of  the  uncoupled 
and  undamped  beam  and  cavity  are  present,  some  variations  occur  in  the  system  of  interest 
since  it  involves  not  only  coupling  between  the  cavity  and  beam  equations  but  also  includes 
damping  in  the  beam. 

In  the  rest  of  the  examples,  harmonic  forcing  functions  of  various  frequencies  but  mod¬ 
eling  a  uniform  pressure  excitation  are  applied  to  the  beam,  and  results  demonstrating 
the  reduction  of  interior  cavity  pressure  and  beam  displacement  under  feedback  control  are 
presented.  Although  the  forcing  functions  were  chosen  to  strongly  excite  various  system 
frequencies,  the  results  show  that  input  of  the  optimally  controlling  voltage  u(t)  uniformly 
reduces  both  the  pressure  and  the  beam  displacement  and  maintains  them  at  a  very  low  level 
of  magnitude  throughout  the  time  intervals  of  interest.  Moreover,  as  indicated  by  the  results 
of  Section  4.3,  this  methodology  is  useful  not  only  for  maintaining  a  low  level  of  pressure 
oscillations  throughout  the  time  interval,  but  also  for  reducing  pressure  oscillations  which 
have  been  allowed  to  build  up  in  the  cavity. 

Numerical  results  have  indicated  that  for  a  1-D  uniform  periodic  forcing  function,  the 
best  results  can  be  obtained  with  one  centered  patch  with  the  amount  of  control  increasing 
with  increasing  patch  length.  By  increasing  the  length  of  the  patch,  one  can  also,  in  many 
cases,  decrease  the  amount  of  voltage  needed  to  reduce  interior  pressure  oscillations.  In 
the  examples  that  we  have  considered  however,  the  magnitude  of  the  controlling  voltage 
has  always  remained  in  a  range  which  is  physically  reasonable  to  apply  to  the  piezoceramic 
patches. 
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